Population-only decay map for n-qubit n-partite inseparability detection 
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We introduce a new positive linear map for a single qubit. This map is a decay only in populations 
of a single-qubit density operator. It is shown that an n-fold product of this map may be used for a 
detection of n-partite inseparability of an n-qubit density operator (i.e., detection of impossibility 
of representing a density operator in the form of a convex combination of products of density 
operators of individual qubits) . This product map is also investigated in relation to a variant of the 
entanglement detection method mentioned by Laskowski and Zukowski. 
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Detection of entanglement in a mixed state is of pub- 
lic interest. There have been many methods proposed, 
such as the partial transpose criterion [l[, the range cri- 
terion 0, the reduction criterion [3,0], the majorization 
criterion [f|, the cross norm criterion Q, the matrix re- 
alignment criterion Q, the entanglement witness crite- 
rion [g,lS|) etc. (There are many review articles on these 
subjects, e.g., Refs. [IS El, 03.) 

Entanglement of a density operator p is defined as in- 
separability of p. For a bipartite system consisting of 
subsystems A and B, p is called separable iff there is a 
decomposition [13( 
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using non-negative probabilities pi (X^Pi = 1). Here, 
variable i is a label. When this decomposition is impos- 
sible, p is called inseparable. 

A multiseparable [Tij state (or fully-separable [To] ] 
state) of n qubits is represented by a density operator 
that can be decomposed in the form of 



(2) 



using non- negative probabilities pi (J2iPi = -0- Any 
density operator that is not multiseparable is called an 
n-partite inseparable density operator of n qubits or an 
n-inseparable density operator of n qubits. We use the 
term n-qubit n-partite inseparability in this meaning. 

To investigate the inseparability, positive but not com- 
pletely positive (PnCP) linear maps are useful pj[ M. [l6| . 
A positive linear map A is a map such that a > ==>■ 
A(er) > for a density operator a (a > means 
(w|o"|i)) > for V|v) in the Hilbert space TL.) For a 
bipartite system whose Hilbert space is Ha ® 7~Cb, if 
we have a density operator p in the form of Eq. |T]), 



iPip\ A ^ <8 A(p[ S ') is kept positive. 
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(iM®AW)(p) = 
By the contraposition of this relation, if (jM ® A^)(p) 
is non-positive, then p is inseparable [f|. Of course, this 
is true also when the map (A^ ® /l B l) is used instead of 
(I^®A^). 

A similar argument is possible for an n-qubit n-partite 
density operator. When a density operator p is in the 
form of Eq. ^j, any positive linear map A acting on a 
single qubit preserves the positivity of (/ <g> A)(p). Here, 
I is the identity map acting on the rest of the qubits. 
When (1(E) A) (p) is a non-positive operator, then p is not 
multiseparable . 

Among such PnCP maps, there are two well-known 
maps: one is a transpose At of a matrix, [Ax(cr)] jj = Oji 
[l[; the other one is the map [|[ Ah (a) = ITrcr — a. 
These maps have been used for the study of entangle- 
ment in bipartite systems intensively. PnCP maps have 
been investigated also for multipartite systems, e.g., the 
depolarizing map [lTj ■ 

Related to the general multiseparability criterion, 
recently, one important theorem was mentioned by 
Laskowski and Zukowski [18j |: 

The Laskowski- Zukowski theorem. Given a bi- 
partite k-separable density operator, the absolute value of 
an antidiagonal element is less than or equal to (l/2) fc . 

As they claimed that this is obvious, the proof is very 
simple: 

Proof. Let us write the Hilbert space of a fc-partite sys- 
tem as using those for subsystems, 
7^[ p j] (j — l ? . . . ; fc). Consider a fc-separable density op- 
erator pIop'- ,Pfel = EiP^f' 1 ® ••• pf k] with Pl > 
such that J2iPi = 1- Any density operator has antidi- 
agonal elements whose absolute values are less than or 
equal to 1/2. Hence, for Vi, pf 1 ^ ® ■ ■ ■ ® pf^ has an- 
tidiagonal elements less than or equal to (l/2) fc . Thus 

picp' "' P ' i ' has antidiagonal elements less than or equal to 
(1/2)*. □ 

By the contraposition of the relation in this theorem, 
we have the corollary: 
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The Laskowski-Zukowski corollary. Given a k- 
partite density operator, if it has an antidiagonal element 
whose absolute value is larger than (l/2) k , the density 
operator is not k-separable. 

This corollary is useful especially in detecting entan- 
glement of n qubits under the condition that all possible 
splittings are considered. When we find an antidiagonal 
element whose absolute value is larger than (1/2)", we 
immediately finish this task. 

In this report, detection of n-qubit n-partite insepa- 
rability is concerned. We first introduce a theorem (a 
variant of the Laskowski-Zukowski corollary) to detect 
n-qubit n-partite inseparability from general off-diagonal 
elements, as a subtopic. Second, we will see the main 
topic where a new positive linear map Ap is defined as a 
single qubit operation. We find that the n-fold product 
Ap n may be used to detect n-qubit n-partite insepara- 
bility. This product map works, under some condition, 
for the class of entanglement that can be detected by the 
new theorem. In addition to these topics, we compare the 
partial- Ap map to the partial transpose and show a few 
examples of detecting n-qubit n-partite inseparability. 

Theorem 1. Suppose that we have an n-qubit n-partite 
density operator p = pi 1 '---*™]. Then, p is not multisepa- 
rable if it has an (a,b) off-diagonal element c a b such that 
\c a b\ > l/2 h ( a,b > ) where h(a,b) is the Hamming distance 
between a and b when expressed in binary notation. 

Proof. Let us write the Hilbert space of an n-qubit system 
as H [1] <g> • • • <g> H [n] by using those for individual qubits, 
W^'l (j = 1, . . . , n). Consider an n-qubit n-separable den- 
sity operator plop"'" 1 = J2i PiPi ® ' ' * ® Z 9 !™' where p t > 
and YliPi — 1- Any density operator of a single qubit 
has diagonal elements [the (0,0) element and the (1,1) 
element] less than or equal to 1 and the antidiagonal el- 
ements [the (0, 1) element and the (1, 0) element] whose 
absolute values are less than or equal to 1/2. Hence for 

Vi, p- 1 ^ ® • • • ® p[™' has (a, b) off-diagonal elements c a b 
for which \c ab \ < (l/2) h( - a ^ is satisfied. This is easy to 
notice when we recall that h{a, b) is the number of dif- 
ferent bits when a and b are compared as binary strings. 

Thus, plop"'"' has (a, b) off-diagonal elements c a b such 
that |c a h| < l/2 h ( a ' b \ By the contraposition of this rela- 
tion, the proof is completed. □ 

Let us look at the illustrating application of the the- 
orem. The theorem is useful especially when we know 
the elements of a density operator. If we find an (a, b) 
element c a ^ whose absolute value > l/2 h ( a,b \ we imme- 
diately find that entanglement is detected according to 
Theorem Q] Indeed, in order to detect that a density 
operator is not multiseparable, using the theorems is of- 
ten more succinct than using ilh—>k-i,k+i,...,n] ^ f or 
several different values of k. Consider the three-qubit 
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(0 < 6 < 1). It is known that this is positive under a par- 
tial transpose with respect to the first qubit (namely, the 
left most qubit when expressed in the binary expression). 
This matrix was originally introduced by Horodecki Q on 
the basis of the range criterion as an inseparable density 
operator with positive partial transposition for a 2 x 4 bi- 
partite system. In the present context, our interest is to 
find it not multiseparable. One of off-diagonal elements 
of p' 1 ' 2 ' 31 is (7|p^ 2 ' 3] |4) = v / T — 62/(146 + 2), which is 

larger than l/2 h ( 7 - 4 ' = 1/4 when 6 < (Vb7 - 7)/4 ~ 

n 2 3i 

0.137. Thus we find that p b ' is not multiseparable 
according to Theorem [1] if 6 < 0.137. In this example, 
the theorem is, of course, weaker than the range criterion 
with which inseparability is detected for < 6 < 1. 

It might be of interest which PnCP map can detect 
entanglement that can be detected by the theorem. In 
the next step, we will see a positive linear map acting 
on a single qubit. According to Propositions 1 and 2 
shown later, the n-fold product of this map works for 
this purpose in either the case where we focus solely on 
antidiagonal elements or the case where all the phase 
factors of individual elements are equal in the left lower 
side (lower than diagonal elements) [or, equivalently, in 
the right upper side (upper than diagonal elements)] of 
a density operator. 

As we have mentioned, we introduce a positive linear 
map acting on a single qubit. 

Definition. The positive linear map Ap for a single 
qubit is defined as follows: 

/ (0H0)+(iHD (0H1) \ 

^ Ap( ^( (iMo) <°h°!Ahi> J ' (4) 

where a is a single-qubit density operator and we use the 
basis {|0>,|1». 

This map may be called a population-only decay map 
because of its definition. 

Note. It is easy to find that Ap is a positive map for a 
single qubit. For a density operator 

\X* Cn J 

the relation |ac| < y/coocn < \ holds. Ap acting on 
a changes diagonal elements into 1/2 and preserves off- 
diagonal elements. Therefore it is a positive map for a 
single qubit. 
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Remark (i). When we see the effect of Ap in the Bloch 
ball of a single qubit, Ap is a projection onto the x-y 
plane as illustrated in Fig. [T] 




FIG. 1: The effect (shown as an arrow) of Ap acting on a 
single-qubit density operator p(x, y, z) = 7/2 + xax + y<JY + 
zgz illustrated in the Bloch ball by its coordinate (x,y,z). 
Here, ax, cry, and oz are the Pauli matrices. The coordinate 
is changed to (x,y, 0) by Ap. 



As is well known, this logic was explained intensively 
by M. Horodecki and P. Horodecki Now we will see 
two remarks relating to this observation. 

Remark (i) . In detecting entanglement of two qubits, the 
partial population-only decay / ® Ap is not so strong as 
the partial transpose 7(8 At, where At is the single-qubit 
transpose operation. For a two-qubit density operator 
/qI 1 ' 2 !, we can write that 

(/^Ap)^ 1 - 2 !) = ±[ p W + (I®A x )(I®A T )(p^)], (6) 

where Ax is a single-qubit NOT operation: A^fcr) = 
XaX (here, a is a 2 x 2 matrix and X = |0)(1| + |1)(0|). 
Thus, (J <£> ApXpI 1 ' 2 !) is sometimes positive for a density 
operator having the property of negative partial transpo- 
sition. It is clear that (I ® AtX/c^ 1,2 ') is non-positive if 
(I (g) Ap)(pl 1 ' 2 l) is non-positive. 

For example, consider the two-qubit isotropic states 



[1.2] 



\M0\ + sl/4 
l + s ' 



Remark (ii) . Let us represent the product of the map Ap 
acting on the fcth qubit and the identity map acting on 
the rest of the qubits by j[i,-.*-l.k+l,...,n] <g This 
product map changes the value of the {i\ . . ■ i n >ji ■ ■ ■ jn) 
element (here, i\, . . . ,i n> j\, . . . , j n € {0, 1}) of an Hermi- 
tian operator £ into: 

*n|C|il ■ • • jfc-lljfc+1 ■ 



(i) UiH 



•ifc-iOifc+i. 
+ {it... i k -ili k+ i . 
when ik=jk, 
(ii) (h . . .i„|£|ji . ..j n ) 
Thus, we have 



• Jn 



■Jn)) 



when i k ± j k . 



,fc-l,fe+i, 



AP)(p) = 



E [ 

a,c,b,d 



-(ac|Tr fc p|bd)(|a0c)(b0d| + |alc)(bld|) 

+ (a0c|p|bld)|a0c)(bld| + (alc|p|b0d)|alc)(b0d|] , 

where Trfc/j is the density operator after tracing out the 
kth qubit; a= i x . ..ik-i, c = i k+1 . ..i n , b = j\ ■ ■ ■ jk-i, 
and d = j k +i ■■■j n - 

Observation. Suppose that we have a multiseparable 
density operator of n qubits, 



[1] 



'Pi 



(5) 



where \<f>) G {^(|00) ± |11)), ^(|01) ± |10>)>; a < -4 or 
< s. It is known that p is ^ is entangled for < s < 2 

[1 21 

because (I (8) At)(p; S q ) has eigenvalues (s — 2)/(4s + 4) 
and (s + 2) /(4s + 4) (with the multiplicity of 3 for the 
latter one) irrespectively of 1 1/>) . When we use the partial 
population-only decay map, we have 



(/0Ap)(pj s 1 ( ; 2 ])e{-± 



4 2(1 + s 
1 ±. 1 



■(|00)(11| + |11>(00|), 
(|01)(10| + |10)(01|)}. 



4 2(1 + s) 

11 21 

Wc find that (7(g) Ap) (p\ m ) has eigenvalues (s — 1)/ (4s + 
4), (s + 3) /(4s + 4), and 1/4 with the multiplicity of 
2 irrespectively of \qb). The first eigenvalue is negative 
for < s < 1. This interval in which entanglement is 
detected is half of that we obtain when I At is used. 

Remark (ii) . It should be remarked that Ap" can be used 
for detecting entanglement. In contrast, if we use the 
single-qubit transpose operation At, A^ n is nothing but 
a transpose and it cannot be used for the detection of 
entanglement. 

For example, consider an obviously entangled pure 
state \ip) = y/p\00) + v / l __ p|H) with < p < I. Then, 



Then, 



.,k-X,k+\ r 



1 ® Ap ' ) (p) is positive. The con- 



traposition of this relation leads to that, for a density op- 



erator p, if (T 



,fc-i,fc+i,...,n] (g)Ap^)(p) is anon-positive 
operator, then p cannot be represented in the form of Eq. 
(0) i.e. p is not multiseparable. Similarly, for a density 
operator p, i/(Ap")(p) is a non-positive operator, then p 
is not multiseparable. 



( 



(I®Ap)(|^|) = 



Vp(i-p) o 






1-P 

2 





Vp(i-p)\ 




1-P / 

2 / 



This has four eigenvalues 
(v/-12p 2 + 12p+l)/4 



p/2,(l - p)/2,l/4 ± 
It has a negative eigenvalue 



4 



when < p < 1. Similarly, for the same pure state, 



These propositions illustrate an ability of . 



it fails 



(A P ® A P )(|¥>)(^|) = 



^/p(l-p)\ 







Wp(i-p) o 



This has four eigenvalues, 1/4,1/4, l/4± y/p — p 2 '. It 

has a negative eigenvalue when \ — VJ < p < 5 + 4 ■ 
Indeed, the entanglement detection using Ap ® Ap docs 
not work for the intervals < p < ^ — ^ — 0.067 and 



VI 
4 



0.933, but still works for the wide 



1 > P > 3 
range of p. 

Using Ap™, we prove the following propositions using 
lemmas shown later. 



Proposition 1. For an n-qubit density operator p — 
p[i,...,n]^ (y\® n )(p) ^ holds if p has an (a,b) antidiago- 
nal element c a b [i.e., h(a,b) = n] such that \c a b\ > 1/2™ 
where h{a, b) is the Hamming distance between a and b 
when expressed in binary notation and a, 6 € {0, . . . , 2™ — 
1}- 

Proof. Consider the n-qubit density operator p. It is as- 
sumed that it has some off-diagonal element c a b such that 
h(a,b) = n and \c a b\ > 1/2". Then apply the operator 
to p. This leads to that 



A ®n 

ilp 



( 7=(Af")(p)=diag(l/2",...,l/2") 

+ off-diagonal elements. 



to map a density matrix to a non-positive Hermitian ma- 
trix for some obviously entangled state, as we have seen, 
but it succeeds for some state positive under a partial 
transpose, such as the state of Eq. (J3j) for t> < 0.137. 

In proving above propositions, we have utilized the fol- 
lowing lemmas. 

Lemma 1. Suppose that an n-qubit density operator p 
has (i,j) elements Cij such that Arg(cy ) = Arg(cfcz) for 
i > j and k > I. Consider the (a,b) element c a b of 
p. Then, (Ap™)(p) has the (a,b) element c a b such that 

\Cab\ > M/2"-^ fc ). 

Proof. Consider binary bits x and y; binary strings a 
and b with bit length k — 1; binary strings c and d 
with bit length n — k. Consider the element of p, 
(accclplbyd). Let p evolve under i[i,—,k-i,k+i,...,n] gjyYp^. 
Then, (axc\(I^>-' k-i,k+i,..., »] ® Ap^X/ZJlbyd) is equal to 
i((a0c|p|b0d) + (alc|p|bld)) if x = y; otherwise it is 
equal to (&xc\p\hyd). Note that if aOc > bOd, then 
ale > bid; if aOc < bOd, then ale < bid. Hence, 
j[i,...,*:-i,t+i,.-,«]g J \f' does not change the argument but 
changes the absolute value of the (axe, hyd) element by 
the factor > 1/2 in the case of x = y, because of equal ar- 
gument values in the (aOc, bOd) and (ale, bid) elements. 
Thus, if we continue to apply T^,-,k-i,k+i,...,n] ^ A W 
from k = 1 to k = n, then for each k, reduction in the 
absolute value by the factor > 1/2 occurs if x — y. There- 
fore \c ab \/2 n - h ^ is a lower bound of \c ab \. □ 



Note that the (a, b) element is unchanged by Ap™ when 
h(a,b) = n. Thus |(a|cr|6)| = |c ah | > 1/2™. Hence a is a 
non-positive operator according to Lemma 2. □ 

Proposition 2. Suppose that we have an n-qubit density 
operator p — p^--> n i such that its off-diagonal elements 
c^ satisfy the relation Arg(cjj) = Arg(cki) for i > j and 
k > I. Then, (Ap")(p) ^ holds if p has an (a,b) off- 
diagonal element c a b such that \c a b\ > l/2 h ( a ' b K 

Proof. Consider an n-qubit density operator p with off- 
diagonal elements Cy such that Arg(cy) = Arg(cfcz) for 
i > j and k > I. It is assumed that it has some off- 
diagonal element c a b such that |c afc | > l/2 h ( a ' b \ Then 
apply the operator Ap™ to p. This leads to 

a =(Af ")(/>) =diag(l/2",... ) 1/2") 

+ off-diagonal elements. 

According to Lemma 1, \(a\a\b)\ > \c ab \/2 n - h< - a ^ > 
1/2™. Hence a is a non-positive operator according to 
Lemma 2. □ 



Lemma 2. Suppose that an Hermitian operator a has 
an (a,b) element s a b G (D (a^b) such that \s a b\ > 1/2™ 
and has diagonal elements su = 1/2™ for V*. Then, a is 
a non-positive operator. 



Proof. Let us use the vector \v) 
Then we have (v|cr|u) = s aa +Sbb- 

\Sab\) <0. 



= |o) - Kb/Kb\)\b). 

\s ab \-\Sab\ = 2(1/2™- 
□ 



To summarize, we have shown a variant of the 
Laskowski-Zukowski corollary that can be used to de- 
tect the n-qubit n-partite inseparability from any one of 
off-diagonal elements of a density operator. To detect 
the n-qubit n-partite inseparability, we have also intro- 
duced and investigated a new positive map Ap acting on 
a single qubit and its n-fold product, Ap™. 

The authors are thankful to Mikio Nakahara for helpful 
discussions and reading through the manuscript carefully. 
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